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Abstract 

In a two orbital double exchange model we formulate Jahn- Teller coupled charge 
density wave in one electron per lattice site limit. Softening of Jahn- Teller phonons 
corresponding to distortion modes Q2 or associated with perfect nesting of Fermi 
surface leads to this instability at low temperature. The gap equation for charge 
density wave state and its dependences on electron-lattice coupling are calculated 
explicitly when any one of the Jahn- Teller modes is excited cooperatively. We find 
that the Q2 distortion mode yields lowest free energy. The calculated heat capacity for 
each mode shows behavior at low temperature. Effect of electron-lattice interaction 
on collective mode, such as amplitude mode, is more pronounced when the excited 
mode is Q2- 
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I. INTRODUCTION 



It is widely accepted that the collective nature of coupling of electrons to the 
underlying lattice degrees of freedom have great impact on extraordinary properties of 
the colossal magneto resistive (CMR) manganites [1]. The parent compound LaMnOs 
contains Jahn- Teller (JT) ions Mn^~^, which is orbitally degenerate. The on-site 
vibrational modes of the MuOq octahedra are breathing mode Qi and the Jahn- 
Teller modes - Q2 (in-plane distortion mode), Q3 (octahedral stretching mode). Due 
to the couphng between the Cg electrons and Q2, Q3 modes the degeneracy of Mn^^ is 
lifted. In the lattice the MnOe octahedra are co-operatively stretched out, resulting 
a C - type orbital and A - type spin antiferromagnet structure [2]. Several recent 
experimental investigations such as high-resolution electron microscopy [3], THz time- 
domain spectroscopy [4] , high- resolution ARPES [5] , non-hnear electrical response [6] , 
x-ray and neutron scattering [7,8] measurements reveal the evidence of charge density 
waves in different composition of CMR manganites. Moreover, ARPES results [5] 
show Fermi surface nesting, which ensures the possibility of density wave formation 
in these materials. It is well known that in a low dimensional one orbital system, 
due to particular geometry of the Fermi surface, electron-phonon interactions lead 
to periodic lattice distortions. This results a ground state characterized by a gap in 
the single particle excitation spectrum at low temperature, which is known as charge 
density wave (CDW) [9]. In CMR manganites the Jahn- Teller vibrational modes 
play major role to shape the properties of these systems compared to the breathing 
mode phonon Qi which docs not alter the symmetry of the MhOq octahedron. The 
distribution of charge density at two Cg orbitals, d^^-yi and d^z^-r^, depends on excited 
JT mode. In presence of these two orbitals, the electron can be either in any one 
of the orbitals or in any linear combination of two [10]. It appears that the perfect 
nesting of the Fermi surface associated with Jahn- Teller distortions may give rise to 
possibilities of occurrence of charge density waves in these systems. Therefore, the 
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resulting charge density wave will be coupled to the orbital order of the system. We 
denote this density wave as Jahn-Teller coupled charge density wave (JTCCDW). 

In the present work, wc investigate the formation of Jahn-Teller coupled charge 
density wave instability in two orbital, double exchange model with one electron per 
site. The paper is organized as follows. In section II we formulate the Hamiltonian for 
Jahn-Teller coupled charge density wave in two dimension and calculate the JTCCDW 
order parameter for different JT vibrational modes. To determine the dominant JT 
distortion mode, we calculate the free energy and hence heat capacity in the JTCCDW 
state. In the same spirit of CDW system, the spectral density functions of the ordered 
ground state are calculated from the amplitude mode response function of JTCCDW 
state generated due to each JT active mode. The results of our calculations are 
presented in section III. Lastly, in section IV we make some concluding remarks. 

II. FORMALISM 

We consider two orbital double exchange model in two dimension, and large 
Hund's rule coupling limit in which the ground state configuration corresponds to 
the itinerant spin being parallel to the core spin at each site. If one compares our 
model with the manganite (e.g. LaMnOs) system the core spins and the itinerant 
electrons may be identified with the t^g (localized) electrons and the (mobile) elec- 
trons of Mn^~^ ion respectively. The kinetic part of the Hamiltonian in momentum 
representation is given by 

^i = E4m, (1) 

k 

where bI = {d\ ,,, d\ f,) with di and ^2 being the Cg - electron annihilation operators in 
two orthonormal orbitals d,j.2_y2 and d^z'^_^.2 respectively. The elements of the hopping 
matrix T are given by Ti_i = —lM[cos,kx + cos/c^], T2,2 = —QM[coskx + cos/c^], 
and Ti,2 = 0.5v^f[ coskx — cos ky]. The transfer hopping integral t is modified as 
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i — t cos(|) because of double exchange interaction [11]. If the t2g spins are assumed 
classical then 9 is the relative angle between two neighbouring spins. Thus, in case 
of ferromagnetic order (FM) t = t. We consider that the A-tjpe spin antifcrro order 
has set in. Due to strong Hund's rule coupling the transport is restricted to spin 
polarizedelectrons in two dimensions only. So the spin degrees of freedom may be 
omitted in the Hamiltonian without loss of generality. 

The Jahn- Teller coupling parts of the Hamiltonian for Q2 and Q3 mode can be 
written as [10] 

Hq^^qYI <52i {dU2i + dldii), (2) 

i 

Hqs ^g^Qzi {d\idu - 4iC?2i ) , (3) 

i 

where g is the electron-JT phonon coupling strength. To diagonalize the kinetic 
part of the Hamiltonian appearing in equation (1) we introduce the new fermionic 
operators c\ and c| as 

1 1 ^k 1 

Cjfc = sm — di^k + cos — d2,k, 

2 , . Ok , . 
Cfc = cos — - sm — d2,k: iv 

with tan6'fc = ^os°\ ^+cos°\ ^) ^ • '^^^ eigen values of the kinetic energy are given by 
An,fe = — cos kx — cos ky — (— l)"'y^cos^ k^ + cos^ ky — cos kx cos ky where n— 1,2. The 
Fermi sea corresponding to Ai^fc is the intersection of the region — 7r/2 < < 7i/2 
with the region — 7r/2 < ky < t:/2 and corresponding to X2,k is the union of the region 
— 7r/2 < kx < 7r/2 (with all values of ky allowed) to the region — 7r/2 < ky < 7r/2 
(with all values of kx allowed). It is evident from the symmetry of the bands that 
at half-filling X^^^ % ~ ~ ^^2(1) {k+Q) ^'-'^ ^ ~ which shows a perfect nesting at 

— * 

wave vector Q. 

The lattice instability in low dimensional system is usually driven by the softening 
of the wave vector Q due to electron-phonon (e-ph) interactions, which is called Kohn 
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anomaly. This is a precursor to the Peierls transition. Peierls instability brings 
about the formation of the CDW with an accompanying periodic lattice distortion 
which in turn gives rise to a band gap at the Fermi level driving the system from a 
conductor to a insulator. In order to find how the JT phonon Q2 or softens to set 
similar instability, it is needed to calculate the phonon self-energy arising from e-ph 
interaction for Q2 or Q3 mode. The phonon response function is determined by the 
phonon Green's function given by [12], 

- r , , • (5) 

TT yu^ -ojg- (^qi^{q, q , (^) 

Here n(g, lu) is the self-energy acquired by the q -th phonon and cUg is the bare phonon 
frequency. The phonon self-energy is given by 

n(g,a;)=47rVx(?,^), (6) 

where x(?;<^) is the density response function. We denote x(g,(^) as X2(?, and 
X3(g, tu) respectively for Q2 and Q3 distorted modes. 

The phonon response function for Q2 and Qs mode can be written as [13], 

k 

X3{q,i^) = H (( P3{q); psi-q) )) a; (7) 

k 

where p2{q) = (4,fc+gC^2fc + 4,fe+gC?ifc) and p3{q) = (rfl.fc+gC^ifc - 4,fc+gC?2fe)- The rcnor- 
malized phonon frequency with the wave vector q — \ Q \ corresponding to the Q2 or 
Qs mode can be determined from the pole of the phonon Green's function, with the 
static response functions X2{Q,0) or X3(Q,0) respectively. It can be shown [14] that 
X2{Q, 0) and xsiQ, 0) diverges at T ^ 0. 

— * — * 

The response function X2{Q, 0) diverges faster than XsiQ, 0) [14]. The divergence 
of the response functions with decreasing temperature indicate the softening of Q2 
and vibrational modes and hence a possibility of Fermi surface instability. The 
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temperature at which the soft mode frequency for Q2 or vanishes is the Jahn- 
Teller coupled charge density wave transition temperature Tcdw corresponding to 
that mode. 

In this report, we formulate the Jahn- Teller coupled charge density wave when 
either Q2 or mode gets excited cooperatively. For this purpose we use (4) and 
second quantized form of the phonon modes given by (fc^*"'^' +h^l^'^) for the wave vector 
q — Q va. equations (2) and (3). By using canonical transformation [15], we derive 
the effective electron-electron interaction mediated by e-ph interaction as, 

^Q2^ ~ ~ X] I 3kk' I sin Ok' (cJ+qC^/ + cJ4/+Q + cJcI,_^q + cJ^qC^,) 

kk' 

X sin 9k {Ck+qcl + c^^c^+q + Ck^l+Q + (^k+Q^k)^ (8) 
Hqs^ — ~ 5Z I 9kk' I COS 9k' {ck'+qcli + cJcI,_^_q + cJcI,^q + cJ^qcI,) 

kk' 

X COS 9k {cIIqcI + c^cl^Q + c^cl^Q + c^+qcI) , (9) 

where gkk' is considered as average interaction strength and gkk' — — \ Qkw |- It is 
to be mentioned here, while deriving the equations (8 ), (9) we have considered only 
terms like c^^^^^c^*^^^ in the electron-phonon interaction. 

Now we write the total Hamiltonian corresponding to Q2 and as, 

^Q2(Q3) = H^j^ Hf^i^'^ . (10) 

In the mean-field approximation, the JTCCDW order parameter (A^) for Q2 and Q3 
modes are respectively defined as, 

Afc = A2(- sin 9k) ^ - Qkk' \ sin 9k' sin 9k 

k' 

{^k'+Q^\' + ^k'^\'^Q + ^k'^\'+Q + ^k'^Q^k')^ (H) 

Afc = A3 cos = - ^ I gkk' I cos 9k' cos 9k 

k' 

i^k'+Q^l' + ^k''^\'+Q + ^k'^lk'+Q + <^fc'+QCfe/), (12) 
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where A2(A3) is the amphtude of the JTCCDW gap generated due to Q2{Q3) mode. 
Therefore, the JTCCDW Hamiltonian can be obtained from (10) as 



ttQ2{Q3) _ V- 3, _ / It U _ 2t 2 N , \ , ^ 2t 2 _ It A \ 

JTCCDW — ^l.fevCfc Cf,_^QCk+Q) -+- 2^ M,k{^k ^k+Q^k+Q) 

k k 

+ ^^k{Ck+Q^\ + ^k^k+Q + ^k^\+Q + ^k+Q^k)^ 



(13) 



while Afc is different for Q2 and mode according to equations (11) and (12). 
To diagonahze the Hamiltonian H^tcc'dw following two simultaneous 

canonical transformations : 



T*fe 1 Vk 2 

Q;i,fe = sin y + cos — c, 



2 ^fc+Q' 



0fc 1 



fe 2 



= COS — Cfc - sm — C;t+Q, 



(14) 



where cos(^) = -i= ( 1 - 



and 



. <^fe 2 , 1 
Q;2,fc = sm y + cos — Cjfc+Q, 

/32,fc = cos y - sm — c^+g, 



(15) 



. These trans- 



wUh co=(f ) = ^ (^1 - and ) = ^ (1 + 

formations yield the quasi particle energies as E"^^ — —Ey} — ^ + 
E^'^ — —E^ — ^J{\\k + A^), which correspond to resulting four bands cti, a2 
and /?2 respectively. The diagonalized Hamiltonian in presence of either Q2 or 
mode can now be written as, 



H 



JTCCDW 



(16) 



The gap equations can be written for Q2 mode as 



and for the mode as 



tanh (/5 /2) tanh (/5 /2) 



+ 



tanh (/3 E'^'12) tanh (/3 £;fcV2) 



+ 



(17) 



(18) 



where for simplicity we have taken | g^k' \— 9- It is clear that both the amplitudes 
A2 and A3 are symmetric in and ky sector of the Fermi surface, but A^. are highly 

anisotropic in /c-space. The effective temperature, T = ksTt enters in (17) and (18) 
through P — l/{kBTt). The free energy of the system can then be calculated as 

^^-^EMl + e-^"^^), (19) 

P k,i 

where El = E^', E^\ E^^ Hence the heat capacity can also be calculated from 
the relation, 

Cv^-T^. (20) 

We like to point out here that the four band Hamiltonian, similar to equation 
(16), is also obtained by Jackeli et al [16] considering only nearest neighbor Coulomb 
correlation in a double degenerate system. To explain the nontrivial observed order 
in LaMnO^, several studies have been reported based on either purely Coulombic in- 
teraction or Jahn- Teller interaction [1] . In the framework of two orbital model, purely 
Coulombic approaches produce correct spin ordered (A-typc spin antiferromagnetic) 
state [17,18], while giving conflicting results regarding the orbital order that co-exists 
with the A-type spin state. According to Hotta et al [19], Jahn- Teller based calcu- 
lations lead to experimentally observed A-type spin and C-type orbital order in a 
model for undoped manganites. Therefore, in the present calculation focussing on 
undoped manganites, we formulate our model based on the Jahn-Teller coupling. 

The robustness of the ordered ground state is usually tested by looking at the 
response of the system to small fluctuations. In case of CDW state, it is well known 
that fluctuations of the phase and the amplitude of the order parameter results in 
the appearance of the collective modes of the system [9] . The amplitude mode is a 
measure of the CDW order parameter. In the present work, we want to calculate 
the spectral density function for the JTCCDW state, from the amplitude response 
function X2,3ilj ^r the Q2 or Q3 modes. The amplitude response function is given 
by, 
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k,k' 



i=l,2 i=l,2 



where 



6**;') = sin6'fc sin6'fc'; /or Q2 - mode 

— cos Ok cos 9k'] for Qs — mode, (22) 

and (^2) fe = (Cfe''^^ ^ c^+Q similar to Nambu operators, fi is the Pauh matrix. 

Within the random phase approximation (RPA), we obtain, 

X2,3 RPAiQ,^) = T. J „ , ^.1 - (23) 



1 -^X2,3(?,'^) 

The spectral density function is related to imaginary part of Xrpa('^), which is given 
by, 

S (0, = Im xh = (l-g;,g3(0^^yVfo;,3(0.^))^ - 

where X2^3(0)'^) ^^^d X2,3(0)'^) respectively the real and imaginary part of 
X2,3(0, a)). By using equation (21) the expressions for X2,3(0>'^) obtained as, 

X^(0,a;)=4sin2^fe x{Q,oj) 

X^(0,(:-) = 4cos2^fc X(Q,'^), (25) 



with 



x(Q,c^)=E 

k 



cos^ 4>k — ^— ^ 9 + cos^ 4 — ^— ^ ^ 



(26) 



III. RESULTS AND DISCUSSIONS 

In previous section we have formulated a theory for the charge density wave state 
in a Jahn- Teller active system. Self consistent solution of equation (17) or (18) yields 
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A2 or A3 respectively. All parameters are expressed in units of t however temperature 
T is expressed as T = ksTt. In figure 1, we show the variation of A2 and A3 (in 
units of t) with temperature T in a FM spin state. It is evident that the magnitude 
of JTCCDW gap decreases with increasing temperature. The magnitude as well as 
transition temperature Tcdw for Q3 mode is larger than that of Q2 mode for ^ = 1. 
The magnitude of A3 is not always larger than A2. This is depicted in figure 2, where 
variation of JTCCDW gap with g is shown for a very low temperature. It is shown in 
figure 2 that Q2 distortion mode opens up JTCCDW gap at a lower value of g than 
that due to Q3 mode and A3 is less than A2 till g < 0.6. With increasing value of g 
both A2 and A3 increases and for values of g greater than 0.6, A2 is less than A3. So 
it is clear that the JT modes have different electron-phonon interaction dependences. 
To determine the dominant JT mode we show the variation of free energy F in FM 
state, with g in figure 3. At very small values of g, e-ph coupling does not alter the 
energy significantly, hence presence of either of the JT modes yields the same energy. 
For higher values of g the JT coupling becomes pronounced and it is evident from 
figure 3 that Q2 distortion mode acquires lowest free energy for any value of g. 

It is well known that low temperature heat capacity (Cy) of CMR manganites has 
a bearing on many fundamental properties like density of states, lattice contribution 
and spin wave stiffness. Recent experiments on insulating LaMnO^ [20] show that at 
low temperature, Cy/T vs plot gives approximately straight lines. In an attempt 
to compare the lattice contributions of the specific heat in JTCCDW state driven 
by Q2 or Q2 mode, we show a plot of Cy/T vs in figure 4 at low temperature. 
Figure 4 shows that at low temperature the behavior is pronounced for both the 
JT modes but the Debye temperatures seem to be different. 

Very recently in Raman scattering experiment on LaMnO^ [21] by Saitoh et al, the 
orbital excitations were observed. According to Brink [22] , the elementary excitations 
of LaMnOs show both the orbital and phonon character due to the mixing of orbital 
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and phonon modes. In the present system, JTCCDW state occurs due to softening of 
either Q2 mode or mode. Here charge densities at different Cg orbitals depend on 
whether Q2 or phonon couphng is active. The mode couples with the electron 
density difference oidx-^-y-^ and d^z'^^r^ orbitals whereas the Q2 mode interacts with the 
charge density difference between bonding and anti bonding orbitals formed by these 
two orbitals [14]. It is expected that the low lying collective excitations of resulting 
JTCCDW will show different g dependence for Q2 and modes. It is evident from 
equation (24), (25) and (26) that spectral density function 5'(0, to) for amplitude mode 
at some symmetry points are due to any one of the JT modes. As example, all spectral 
density at symmetry point (0, 0) is due to mode while at (tt, 0) only Q2 mode 
contributes to 5'(0, a;). It is also obvious that at these symmetry points the CDW 
gap exists due to the JT mode which contributes to S'(0,d;). However, the 5'(0, cu) as 
well as Afe at (tt, 7r/2) is finite for both Q2 and mode. In figure 5 and figure 6 we 
show 5'(0, a;) at (tt, 7r/2) for different values of g in Q2 and excited CDW state 
respectively. With increasing g the peak shifts to higher values of frequency ui for 
both the JT modes, because the amplitude of the JTCCDW gap increases for both 
the modes. It is important to note that the peak position of 5'(0,a;) shifts more for 
Q2 mode than that of Q3 with equal increase of g. So the electron lattice interaction 
has more pronounced effect on the collective excitations of Q2 active CDW state than 
that of mode. 



IV. CONCLUSIONS 

The charge density wave may have its origin due to softening of the breathing 
mode phonon or Jahn-Teller phonons. Emergence of the Jahn-Teller coupled charge 
density wave is relevant in a system where Jahn-Teller phonons play a crucial role to 
determine the orbital ordering. In the limit of one electron per site and strong Hund's 
rule coupling, we provide a comprehensive study of the formation of charge density 
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wave emphasizing on Cg electron - Jahn- Teller coupling in a orbital ordered system 
with two kinds of orbitals alternating on adjacent sites in the x—y plane. The in-plane 
distortion mode - Q2 is the dominant mode and it opens up JTCCDW gap at smaller 
electron-lattice interaction strength than the octahedral stretching mode - Q3. The 
heat capacity shows behavior at low temperature. At some of the /c-points, the 
effect of the electron-lattice interaction on collective mode such as amplitude mode 
is more pronounced when the excited Jahn- Teller mode is Q2- 

The present calculation is done for a system with one electron per lattice site and 
either of the JT modes is active. The divergence of the phonon response functions 
corresponding to Q2 or mode ensures their softening and hence the formation of 
CDW state. However, in presence of both the modes, as in real systems, it is compli- 
cated to formulate the charge density wave instability from first principle calculation. 
With hole doping in the present system there will be decrease of Jahn-Teller active 
Mn^+ ions and the perfect nesting of the Fermi surface will be destroyed. These are 
expected to reduce the effective electron- JT interaction as well as JTCCDW gap A^, 
which may yield a orbitally disordered system. The work along these lines are in 
progress. 
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Figure captions : 

Fig. 1. Thermal variations of the A2 (sohd hne) and A3 (dashed hne) in FM state, for ^ = 1.0. 
Here, f = ksTt. 

Fig. 2. Variations of A2 (sohd hne) and A3 (dashed hne) with g at f ^ {kBTt) = 10"^ in FM 
state. 

Fig. 3. Variation of free energy F in FM state for Q2 distorted CDW state (sohd hne) and 
Qs distorted CDW state (dashed hne) with ^ at T = 0.005. 

Fig. 4. Cy/T vs plot in FM state for Q2 active (sohd hne) and Q2, active (dashed hne) 
CDW state with g = 1.0. Here Cy is in arbitrary units and T = kBTt. 

Fig. 5. Spectral density functions 5(0,0;), for the JTCCDW amplitude mode as a function of 
frequency uj (in unit of t) in presence of Q2 distortion, at /c-point (tt, 7r/2), for different 
values oig^ 0.45, 1.0, 1.5. 

Fig. 6. Spectral density functions 5'(0,a;), for the JTCCDW amplitude mode as a function of 
frequency ui (in unit of i) in presence of Q3 distortion, at /c-point (tt, 7r/2), for different 
values of ^ = 0.45, 1.0, 1.5. 
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